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Density functional calculations of electronic structure, total energy, structural distortions, and 
magnetism for hydrogenated single-layer, bilayer, and multi-layer graphene are performed. It is 
found that hydrogen-induced magnetism can survives only at very low concentrations of hydrogen 
(single-atom regime) whereas hydrogen pairs with optimized structure are usually nonmagnetic. 
Chemisorption energy as a function of hydrogen concentration is calculated, as well as energy barriers 
for hydrogen binding and release. The results confirm that graphene can be perspective material 
for hydrogen storage. Difference between hydrogenation of graphene, nanotubes, and bulk graphite 
is discussed. 

PACS numbers: 73.20.Hb, 71.15.Nc, 81.05.Uw 
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I. INTRODUCTION 

Discovery of graphene, the first truly two-dimensional 
crystal, and its exotic electronic properties (for review, 
see Refs. Il|l2||3j) initiates a huge growth of interest to car- 
bon materials. Most of activity is focused on electronic 
transport phenomena in graphene, keeping in mind po- 
tential applications for carbon-based electronics. How- 
ever, chemical physics of graphene is also very interest- 
ing, in particular, due to opportunity to use graphene 
for chemical sensors with extraordinary sensitivity^. An- 
other interesting direction of investigations is a possible 
use of graphene for hydrogen storage. One could expect 
that two-dimensional systems could be very convenient 
for this aim. 

In general, carbon-based systems are among the most 
attractive objects for hydrogen storage^. A promis- 
ing storage properties of single-wall carbon nanotubes 
(SWCNT) were first reported in Ref. 0. In last few 
year graphene was used as a model system to study 
the electronic structure and adsorption properties of the 
SWCNT— After discovery of real graphene several 
works appeared theoretically studying the hydrogen ad- 
sorption on graphene, as a special material (see, e.g., 
Refs. l9llIof). It is commonly accepted nowi£&i£ that 
the chemisorption of single hydrogen atom on graphene 
leads to appearance of magnetic moments in the system. 
The magnetic interactions between the hydrogen atoms 
placed at large distances on graphene have been calcu- 
lated in Ref. |9|. However, energetics of various hydrogen 
configurations taking into account carbon sheet relax- 
ation was not studied yet in detail. In earlier works, 
only a very special structure of hydrogenated graphene, 
with all hydrogens sitting on the same side was discussed. 
Here we will demonstrate that actual energetically favor- 



able structure with hydrogenization of the both sides has 
quite different properties and, in particular, turns out to 
be nonmagnetic. 

Earlier a similar structure has been discussed for the 
case of SWCN T 11 : 12 . However, in contrast with the 
SWCNT in graphene there is no specific potential bar- 
rier for hydrogen atoms^ 2 - since both sides of graphene 
are equally achievable for the adsorption which makes 
the situation different. Deeper understanding of the case 
of graphene will be useful also to discuss hydrogen stor- 
age capacity of nanotubes 7 ^ or nonporous carbon 1 ^, as 
well as corresponding experimental results for graphite 14 . 
Effect of curvature on the hydrogen chemisorption in 
fullerenes and nanotubes has been considered earlier in 
Ref. [H. 



II. CHEMISORPTION OF SINGLE HYDROGEN 
ATOM 

To model the hydrogen chemisorption we use a peri- 
odic supercell of graphene containing 32 carbon atoms 
per each hydrogen atom, similar to Ref. 7J. To con- 
sider hydrogen pairs, we will use supercells with 50 car- 
bon atoms for close pairs (neighboring positions of hy- 
drogen) and 72 carbon atoms, otherwise. The density- 
functional theory calculations were performed using the 
SIESTA cod e 16 : 17 which was successfully applied before 
to describe hydrogen on graphene^. We used the same 
technical parameters of the calculation as in Ref. ©. 

To discuss chemisorption on graphene it is worth to 
remind its basic electronic structure. Originally, carbon 
has two 2s and two 2p electrons. These four electrons 
produce different kinds of sp-hybridized orbitals 1 ^. In 
graphene every carbon atom is bounded with three other 
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FIG. 1: (color online) Band structure of a single graphene 
layer. Solid red lines are a bands and dotted blue lines are tt 
bands. 





FIG. 2: (color online) Crystallographic structure of graphene. 
Red and blue circles show A and B sublattices, respectively. 
Labels show the distance from Ao carbon atom (coordination 
sphere numbers). All bonds in graphene are equivalent, the 
double bonds are marked for convenience of comparison with 
other pictures. 



carbon atoms via sp 2 hybridization. There are three a 
orbitals placed in the graphene plane with angle 120° 
and one it orbital along Z axis in perpendicular direction. 
Figure [1] shows the band structure of pure graphene, with 
three a bands laying about 3 eV above and below the 
Fermi level, and 7r band. In diamond, all carbon atoms 
are connected via sp 3 hybridization with four a bands 
separated by a big gap. Breaking 7r bonds and producing 
additional a bond and, thus, transition from sp 2 to sp 3 
hybridization is the main mechanism of chemisorption 
on graphene. The crystallographic structure of graphene 
with two sublattices is shown in Fig. O In pure graphene 
the sublattices are equivalent, but if we bind one of car- 
bon atoms (for example, A in Figured]) with hydrogen 
we automatically break this equivalence. 

To check the computational procedure, we repro- 




FIG. 3: (color online) Picture of local distortions of graphene 
at chemisorption of: (a) single hydrogen atom (Ao); (b) two 
hydrogen atoms bonded with carbon atoms from the same 
sublattice (A0-A2); (c) two hydrogen atoms bonded with 
neighboring carbon atoms from the same side of graphene 
sheet (A0-B1); (d) two hydrogen atoms bonded with neighbor- 
ing carbon atoms from both sides of graphene sheet(Ao-B'i). 
Red and blue circles are carbon atoms from two sublattices, 
white circles are hydrogen atoms. 



duce first known results ^ 8 ' 9 ' 10 for single hydrogen atom 
chemisorbed on graphene. In agreement with the previ- 
ous calculations we have found hydrogen-carbon distance 
about 1.1 A, and shift of the carbon atom bonded with 
the hydrogen one about 0.3A along Z direction. One 
should stress, additionally to the previous results, that 
the atomic distortions are not negligible also for the sec- 
ond and third neighbors of the hydrogen-bonded carbon 
atom Ao (see Figured). Amplitude of the modulation 
of graphene sheet in the perpendicular direction around 
the hydrogen atom was estimated as 0.4A, which is com- 
parable with the height of intrinsic ripples on graphene of 
order of 0.7A found in atomistic simulations^. The ra- 
dius of the distorted region around hydrogen atom turned 
out to be about 3.8A. 

Transformation of the sp 2 hybridization of carbon in 
ideal graphene to the sp 3 hybridization in hydrogenated 
graphene results in a change of the bond lengths and an- 
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TABLE I: 



Dependence of magnetic moments M (in [Mb), chemisorption energies E c h em (in eV), and geometrical parameters (see Figure 
[3]), in degrees and A, on configuration of hydrogen (see Figure [3J; d are interatomic distances and h are heights of atoms from 
graphene plane. 



Configuration 


M 


Echem 


h-AO 


h B i 


h,A2 


angle(C-C-H) 


angle(C-C-C) 


dc-H 


dc-c 


A 


1.0 


1.441 


0.257 


-0.047 


-0.036 


101.3 


115.4 


1.22 


1.496 


Ao-A 2 


2.0 


1.406 


0.285 


-0.040 


-0.096 


102.7 


116.6 


1.132 


1.483 


A0-B1 


0.0 


0.909 


0.364 


-0.088 


-0.069 


102.2 


117.5 


1.077 


1.491 


Ao-B'i 


0.0 


0.540 


0.298 


-0.027 


-0.035 


105.1 


106.7 


1.112 


1.512 



gles. A typical bond length for sp 2 C-C bonds is 1.42A for 
graphene and graphite and 1.47A for other compounds, 
and the standard bond angle is 120°. For sp 3 hybridiza- 
tion, the standard value of C-C bond length is 1.54A, 
and the angle is 109.5°. A typical value for the single 
C-H bond length is 1.086A. One can see in Table U that 
for single hydrogen atom the C-H bond length is close to 
the standard value, but C-C-H and C-C-C angles are in- 
termediate between 90° and 109.5° and 120° and 109.5°, 
respectively. Also, the length of C-C bond is in between 
1.42 A and 1.54A. This means an intermediate character 
of the hybridization between sp 2 and sp 3 . 

A pictorial view of the reconstruction of chemical 
bonds, with the breaking of double C=C bond and forma- 
tion of single C-H bond, is shown in Fig. 0J For the case 
of single hydrogen atom (Fig. this releases two un- 
paired electrons. One of the electrons forms a new bond 
with hydrogen whereas the other is unpaired. The latter 
is delocalized in some rather broad area on lattice^. As 
a result, carbon becomes magnetic (see the Table and 
hydrogen atom also possess a small magnetic moment 
about 0.12 u,B- In general, at the chemisorption of single 
carbon atom the hybridization is still rather close to sp 2 . 
One has to consider another opportunities which can lead 
to sp 3 bonding and possible gain in the chemisorption en- 
ergy. 



III. HYDROGEN PAIRS ON SINGLE-LAYER 
GRAPHENE 

There are four kinds of hydrogen pairs on graphene: 
hydrogen atoms can be bonded by carbon atoms from the 
same sublattice of from different sublattices, on one side 
from the graphene sheet, or from both sides. We use the 
primed indices for the later case. Computational results 
for chemisorption of hydrogen pairs are presented in Fig. 
[3] and in Table |TJ Chemisorption energy per hydrogen 
atom for the case A0-A2 (next-nearest-neighboring car- 
bon atoms, both hydrogen atoms are from the same side) 
is not significantly different from that for single hydro- 
gen, whereas chemisorption by carbon atoms from differ- 
ent sublattices turns out to be much more energetically 
favorable. 

To understand the difference, one has to study what 
happens with chemical bonds in all these cases. In Figure 




FIG. 4: (color online) Sketch of chemical bonds for chemisorp- 
tion of hydrogen: (a) single hydrogen atom (Ao); (b) two 
hydrogen atoms bonded with carbon atoms from the same 
sublattice (A0-A2); (c) two hydrogen atoms bonded by near- 
est carbon atoms (A0-B1); (d) two hydrogen atoms bonded by 
next-nearest carbon atoms from different sublattices (A0-B3). 
Big red (dark) and blue (light) circles are carbon atoms from 
different sublattices, small white circles are hydrogen atoms, 
small black circles are unpaired electrons. 

SJd, we can see that for the case A0-A2 the situation is 
basically the same as for the single hydrogen, namely, two 
broken bonds produce two unpaired electrons with strong 
ferromagnetic coupling between their spins (dependence 
of the exchange interactions from interatomic distance 
was studied in detail in Ref. These electrons in the 
Aq-A u case are not paired and produce new chemical 
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FIG. 5: (color online) Partial densities of states for carbon 
atom bound with hydrogen (solid red line), hydrogen atom 
(dashed green line), and distant carbon atom (from fourth 
coordination sphere) (dotted blue line) for the case Ao-Bi' 
configuration (see Fig. [3J1). 



bonds, the bond distances and angles for A -A 2 being 
intermediate between those typical for the sp 2 and sp 3 
hybridization (see Table UJ. 

The situation Ao-Bi is essentially different. One can 
see from Fig. 0t that, when the double bond between 
Ao and Bi carbon atoms transforms into the single one 
two unpaired electrons appears and both of them partic- 
ipate in formation of covalent bonds with the hydrogen 
atoms. For the case of more distant carbon atoms, say, 
Ao and B3 we can see a similar situation (Fig. QJ1). Cor- 
responding changes in the electronic structure for this 
case is displayed in Fig. The density of states for car- 
bon atoms bonded with hydrogen in redistributed, de- 
creasing in the region between -2.5 to 5 eV (the energy 
is counted from the Fermi level) and increasing near ±7 
eV. These changes correspond to a transition from sp 2 to 
sp 3 hybridization which makes graphene-like electronic 
structure more "diamond-like" transforming the 7r band 
crossing the Fermi level (see Fig. [l| to fourth a band ly- 
ing far from it. At the same time, the electronic structure 
of fourth neighbors are very close to electronic structure 
of pure graphene (Fig. [5]). fn the case of chemisorp- 
tion by carbon atoms from different sublattices there are 
no unpaired electrons and no magnetism. In the work 9 
this situation was described as antiferromagnetic which 
is not quite accurate as we believe. Actually, the local 
magnetic moments just do not survive in this case. The 
absence of unpaired electrons and broken bonds leads to 
chemisorption energy gain in comparison with the A-A 
case described above. 

Additionally, we can see in Table U that the C-C bond 
length for the case A -Bi is close to the standard one for 
sp 3 hybridization. However, the bond angles are closer 
to those for sp 2 hybridization, and the chemisorption en- 
ergy for the case Ao-Bi is higher than for Ao-B'i. To 
understand the difference, one has to investigate struc- 



tural distortions of graphene sheet. Chemisorption of 
hydrogen by Ao carbon atom induces its shift up per- 
pendicular to the plane, together with shifts of atoms 
Bi and A2 in the opposite direction. The chemisorption 
on carbon Bi atom shifts Bi atom up and Ao and A2 
atoms down. Therefore, for the case Ao-Bi both Ao and 
Bi carbon atoms move simultaneously in the same direc- 
tion. As a result, the bond angles become close to those 
typical for sp 3 hybridization. On the contrary, in Ao- 
B'i case the chemisorption of hydrogen from the bottom 
by Bi carbon produces shifts up for Ao and down for Bi 
carbon atoms that coincide with the lattice distortion for 
the bonding of hydrogen by Ao from the top. In the case 
A -B'i the lattice distortions produced by chemisorption 
of each hydrogen atoms are consistently working in the 
same direction providing the lowest chemisorption energy 
and bond lengths and angles closest to the standard ones 
for sp 3 hybridization (see Table UJ). 

The calculated dependence of the chemisorption en- 
ergy on the distance between carbon atoms bonded with 
hydrogen is presented in Fig. O One can see that for 
all types of pairs the chemisorption energy for the hy- 
drogen atoms closer than 5A is lower than for larger 
distances. Independently on the distance, the nonmag- 
netic A-B pairs are more energetically favorable than A-A 
pairs and than noninteracting hydrogen atoms. One can 
assume therefore that observation of hydrogen-induced 
ferromagnetism 9 is possible only for a very low concen- 
tration of hydrogen when the distance between hydro- 
gen atoms are higher than 12A. Our results seem to be 
in a qualitatively agreement with the experimental data 
on hydrogen chemisorption on highly-oriented pyrolitic 
graphite (HOPG)^. The pairs A0-B3 have been observed 
which correspond to minimal energy for the one-side hy- 
drogenation of graphene, according to our results (see 
Fig. [B]). Also, at hydrogenation of fullerenes Ceo the 
pairs Ao-Bi and A0-B3 (1,2 and 1,4, according to chem- 
ical terminology) are usually observed (see, e.g., the re- 
view [24] and references therein) . Instability of magnetic 
state was observed experimentally for C60H24 2 —. Recent 
theoretical results for chemisorption on singl o 21 ' 22 and 
multiple- wall2£ carbon nanotubes are qualitatively simi- 
lar to our results for graphene. 



IV. HYDROGEN CHEMISORPTION ON 
BILAYER GRAPHENE 

Let us consider now hydrogen chemisorption on 
graphene bilayer. We studied the chemisorption of sin- 
gle hydrogen atom and pairs of hydrogen atoms placed 
on one and both sides of the bilayer. The calculations 
have been performed for two different concentration of 
hydrogen, that is, low (32 carbon atoms in each layer 
per hydrogen atom) and high (8 carbon atoms in each 
layer per hydrogen atom). Lattice distortions induced 
by the hydrogen turned out to be different for the case 
of single-layer and bilayer graphene. Whereas the shift 
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FIG. 6: (color online) Energy of hydrogen pair (per atom) 
counted from the energy of single hydrogen atom as a function 
of interatomic distance: Ao-B„ - solid red line with crosses, 
Ao-B' n - dashed green line with crosses, Ao-A n - dotted light 
blue line with filed squares, Ao-A' n - dot-dashed violet line 
with empty squares. 



of carbon atom bound with hydrogen is rather similar in 
both cases, atomic displacements for the neighboring car- 
bon atoms are much smaller in the case of bilayer. This 
is not surprising since interlayer coupling tends to make 
graphene more flat, e.g., sheet corrugations are smaller 
for suspended bilayer membrane than for the single-layer 
one2& 

Computational results are presented in Table [TTJ One 
can see that for low hydrogen concentration the differ- 
ence of chemisorption energies between single hydrogen 
atoms and the pairs is smaller than for the case of single- 
layer. There are two configurations which have a very 
close values of the energy for low concentration of hy- 
drogen, Ao-Bi and A0-B3. For the higher concentration, 
the latter configuration becomes essentially more stable 
since the lattice distortions are more homogeneous in this 
case. The effective interactions between hydrogen atoms 
is more short-range in the case of bilayer and already 
for the configuration A0-B5 the chemisorption energy is 
almost equal to that of two single hydrogen atoms. 

In the case of single-layer the hydrogen positions on 
different sides of the graphene sheet are essentially more 
favorable than those on the same side. Contrary, for the 
case of bilayer this energy difference is small. 



V. HYDROGEN STORAGE PROPERTIES OF 
GRAPHENE 

Chemisorption energy per hydrogen atom for the most 
favorable case of Ao-B'i pairs presented in Table [T] is not 
very high. Another limiting case with much higher ad- 
sorption energy per hydrogen atom corresponds to the 
case of fully hydrogenated graphene which is close to a 
hypothetical compound graphane^I. For the latter case, 



we have found bond lengths 1.526A for C-C bonds and 
l.HOA for C-H bonds, and bond angles 102.8° and 107.5° 
for C-C-C and C-C-H angles respectively, in a good agree- 
ment with the results of Ref. |27|. The calculated values 
are close to the standard ones for sp 3 hybridization, that 
is, 1.54A for the length of C-C bonds and 109.5° for all 
angles. Value of C-H bonds are also very close to the 
standard 1.09A. 

We studied transition from single pairs to complete 
coverage changing the supercell size. The dependence of 
the chemisorption energy on the hydrogen concentration 
is shown in Fig. For fully hydrogenated graphene 
the mass percentage of hydrogen (gravimetric energy 
density), is 7.8 which is over the target value of DOE 
(United States Department of Energy) 6.5 6 . Another 
relevant characteristics for hydrogen storage are energy 
barriers which are necessary to overcome to start hy- 
drogenation and dehydrogenation. They correspond to 
the chemisorption energy per hydrogen atom for singe 
hydrogen pair and for fully hydrogenated graphene, re- 
spectively. We found for these quantities 0.53 eV (25.5 
kJ/mol) and 0.42 eV (20.3 kJ/mol). The latter value is 
close to the experimental one, 19.6 kJ/mol, for hydrog- 
enized nanotubes 6 . These values look quite reasonable 
in view of potential applications of graphene for the hy- 
drogen storage. Transformation of electronic structure 
with increasing hydrogen concentration presented on the 
insets of Fig. [7] Minimal mass hydrogen concentration 
which results in opening of energy gap at the Fermi level 
is about 4.04 (50% coverage), the gap value being 1.75 
eV. This seems to be, potentially, an interesting predic- 
tion for experiment, although it is not clear whether it is 
possible to stabilize this configuration or not. 

The computational results under discussion have been 
obtained in the generalized gradient approximation, 
GGA, which is a common practice for electronic structure 
calculations of H-C systems^^^. To estimate possible 
errors we have calculated the desorption energy in the 
local density approximation (LDA) as well. We have ob- 
tained the value 0.62 eV, in comparison with the GGA 
result 0.42 eV so the difference is essential. In more de- 
tail, the question was studied in Ref. [H with the conclu- 
sion that GGA is more reliable than LDA for this kind 
of problems. 

In previous worksZ^ hydrogenation of graphene was 
studied as a model of that of SWCNT. However, these 
two situations are not identical due to curvature of the 
nanotubes. In Fig. [5]we sketch the SWCNT, value h from 
Fig. [HJd corresponding to the sum of the values Hao and 
h A 2 from Fig. Hi. In SWCNT h = a 2 /2R, where a is the 
lattice parameter for graphene, 2.46A, and R is the radius 
of nanotube. Typical diameters of the SWCNT 10-j-15A 
correspond to the values of h from 0.605 to 0.375A. At the 
same time, single hydrogen atom on graphene produces a 
distortion with the value h = 0.293A, that is lower than 
for the SWCNT of standard diameter. This value is close 
to those for theoretical estimations of maximum of the 
SWCNT diameter, 41.6^ and 49.9^A. 
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TABLE II: 



Chemisorption energy E c h em per hydrogen atom (in eV), height h of carbon atom bound with hydrogen up to the layer, 
and interlayer distance d (in A) for graphene bilayer for different hydrogen concentrations and configurations of chemisorbed 
hydrogen. 



Concentration 


Configuration 


E 1, 


h 


d 


Low 


A 


1.28 


0.639 


3.237 




Ao-Bi one side 


0.715 


0.570 


3.222 




Ao-Bi both sides 


0.713 


0.615 


3.149 




A0-B3 one side 


0.720 


0.477 


3.237 




A0-B3 both sides 


0.733 


0.453 


3.237 


High 


Ao-Bi one side 


0.885 


0.445 


3.174 




Ao-Bi both sides 


0.850 


0.426 


3.041 




A0-B3 one side 


0.381 


0.359 


3.262 




A0-B3 both sides 


0.390 


0.349 


3.198 



0.6 




Hydrogen weight, % 



FIG. 7: (color online) Dependence of the chemisorption en- 
ergy per hydrogen atom on the mass hydrogen concentration 
(gravimetric energy density) . The insets shows total densities 
of states for for (a) 2.06 and (b) 4.04 mass hydrogen concen- 
tration. 



On the other hand, multiple-wall carbon nanotubes 
(MWCNT) have typical diameters about 50A and bilayer 
graphene can be a reasonable model to study hydrogena- 
tion of the MWCNT. Moreover, partial graphitization 
and presence of metallic catalysts strongly influence on 
adsorption properties of SWCNT- whereas graphene is 
perfectly pure material. Another problems for hydro- 
genation of nanotubes are how to provide an access of 
hydrogen to their surface in an array^i and high enough 
flip-into energy barrier—. Carbon Is X-ray photoemis- 
sion spectra (XPS) of the SWCNT before hydrogenation, 
after hydrogenation, and after dehydrogenation reported 
in Ref. [H are all different that could be in part due to 
defect formation whereas graphene has a very high va- 
cancy formation energy (up to 8 eV) which means much 
higher stability of graphene under high temperatures and 
pressures. 

At last, we compared hydrogen storage properties 
of graphene and graphite nanofibers (GNF), that is, 
very small graphite platelets, with a size of order of 
30-^500A22. Raman spectra for graphene multilayers 



a) 




FIG. 8: (color online) (a) Position of carbon hexagons on sur- 
face of SWCNT. Red (dark) and blue (light) circles are carbon 
atoms from different sublattices. (b) Position of carbon atoms 
on radii of SWCNT. 



become very close to ones for graphite when number 
of layers is five or more so one can assume that five- 
layered graphene is already similar to the bulk case^. 
To model the GNF we used therefore five-layer graphene 
slab. Complete one-side hydrogenation of GNF, as well 
as of graphene, is impossible and only 50% hydrogena- 
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TABLE III: 

Dependence of chemisorption energy (in eV), interlayer distance d, and geometrical parameters in A (see Fig. [3} on numbers 
of graphene layers for 50% hydrogenation of one side of the top layer. 



Number of layers 




d 




Kb 


dc-H 


dc-c 


1 


1.775 




0.106 


0.143 


1.158 


1.475 


2 


1.452 


2.88 


0.142 


0.198 


1.154 


1.475 


5 


1.621 


3.124 and 3.353 


0.133 


0.116 


1.164 


1.468 



tion of the top layer is supposed to be the maximum 
(all carbon atoms from one of the sublattices are bonded 
with hydrogens) that corresponds to approximately 1% 
of gravimetric energy density. The calculated chemisorp- 
tion energy per hydrogen for single-layer graphene with 
the same gravimetric energy density is 0.32 eV lower than 
for five-layer graphene. The maximum load of the five- 
layer graphene is 2% of the gravimetric energy density 
that is about four times smaller than for the single-layer 
graphene. Results of calculation for case of 50% hy- 
drogenated surface of graphene single-layer, bilayer and 
graphite (five-layer of graphene) are presented in Ta- 
ble IIIII For all three structures chemisorption energies, 
structural changes, magnetic properties and electronic 
structures are essentially different. Some differences, e.g., 
in the length of C-H bond are negligible, but many others 
(amplitude of bending distortions) are significant. De- 
tailed comparison of chemical and structural properties 
of single-layer and multilayer graphene will be reported 
elsewhere. 



ics of different hydrogenated grcphenc layers. Our results 
support a suggestion that graphene may be a promising 
material for the hydrogen storage. Equivalence of two 
sides of graphene distinguishes it drastically from the 
nanotubes. We have shown that the most stable config- 
uration of low hydrogenated grephene layer corresponds 
to the non-magnetic pair hydrogen atoms attached to 
the different A-B sublattices of graphene from the dif- 
feren sides. It is worth to emphasize that single-layer— 
or bilaye r 35 ' 36 graphene should be rather carefully used as 
models of structural and chemical properties of graphite 
and its derivates. Also, a comparison of experimental 
results for graphite^ with the computational results for 
graphene sheets^ requires additional investigation. 
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